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In response to W. G. Hoover's comment >rXiv:1204.0312y2] on our work arXiv:1203.5968 , we 
show explicitly that the divergence of the velocity field associated with the Nose-Hoover equations is 
nonzero, implying that those equations are not volume preserving, and hence, as often stated in the 
literature, are not Hamiltonian. We further elucidate that the trajectories {q(t)} generated by the 
Nose-Hoover equations are generally not identical to those generated by Dettmann's Hamiltonian. 
Dettmann's Hamiltonian produces the same trajectories as the Nose-Hoover equations only on a 
specific energy shell, but not on the neighboring ones. This fact explains why the Nose-Hoover 
equations are not volume preserving. The Hamiltonian that we put forward with arXiv:1203.5968 
instead produces thermostated dynamics irrespective of the energy value. The main advantage 
of our Hamiltonian thermostat over previous ones is that it contains kinetic energy terms that 
are of standard form with coordinate-independent masses and consequently is readily matched in 
laboratory experiments. 

PACS numbers: 02.70.Ns, 05.40.-a 



Iii his comment [lj to our preprint [i[ W. G. Hoover 
criticizes our statement that the Nose-Hoover equations 
are not Hamiltonian, and maintains that, on the con- 
trary, they are Hamiltonian. We found the arguments 
given by Hoover in support of his statement incorrect. 

Hoover [lj begins with Nose's Hamiltonian of a har- 
monic oscillator of unit mass and unit frequency: 

H Nos6 = (l/2)[(p/s) 2 + q 2 + (p s / T ) 2 ]+Tlns (1) 

where r and T are parameters controlling a time scale 
and the temperature of the thermostat, respectively. The 
corresponding equations of motion: 

q = p/s 2 ; p= -q ; s = p s /r 2 ; p s = p 2 /s 3 - T/s , (2) 

are by definition Hamiltonian Q. Hoover then proceeds 
with a "time-rescaling" to obtain the new set: 



p/s ; p 



-sq 



sp s /r 2 ; p s = (p/s) 2 - T. (3) 



This new set is evidently non-Hamiltonian. A cru- 
cial property of a Hamiltonian flow is that the diver- 
gence of the velocity field in phase space is null (this 
is equivalent to saying that the flow is incompressible, 
i.e., obeys Liouville theorem). Let v be the velocity 
vector field v = (q,p,s,p s ). From follows that 
divv = dq/dq + dp/dp + ds/ds + dp s /dp s = p s /t 2 ^ 0, 
implying that the set of Eqs. © is not Hamiltonian. 

In the next step Hoover makes the change of variables 
v = p/s and £ = p s /r 2 to obtain the well known Nose- 
Hoover equations: 



v : v 



-q-(v ; C=[v 2 -T]/t 2 



(4) 



The divergence of the associated velocity field v = 
(q, v, Q is nonzero: divv = dq/dq + dv/dv + dC,JJ)C, = 
— ( =/= 0, implying that, as we state in our paper [2j, the 
Nose-Hoover Equations (jj) are not Hamiltonian. This 



is indeed a statement that often appears in the pre- 
vious literature. Klages Q repeatedly states that the 
Nose-Hoover equations are not Hamiltonian. Kusnczov, 
Bulgac and Bauer in Ref. Q comment about the Nose- 
Hoover Equations that: " These equations of motion no 
longer retain a Hamiltonian structure. In the extended 
phase space the variables q, p and £ are no longer canon- 
ical and the symplectic structure is lost" . 

Hoover [l| continues with deriving the Nose-Hoover 
equations ((4]) from the Hamiltonian of Dettmann Q : 

ffDcttmann = p 2 /2s + sq 2 /2 + s P 2 /2t 2 + sT In s (5) 

The corresponding canonical equations 

q = p/s 

p = -sq 

s = sp s /T 2 
p s =p 2 /2s 2 ~ q 2 /2-p 2 J2t 2 -T\ns-T 

(6) 

are, by definition, Hamiltonian. Hoover [lj then chooses 
a special value of the Hamiltonian, -ffDcttmann = 0, and 
uses this numerical value of the Hamiltonian in (|6]l. to 
obtain the set of Eqs. <j3j> - We note that replacing the 
Hamiltonian function by the numerical value of the en- 
ergy, as it is determined by the initial conditions, actually 
leads to equations of motion which no longer obey Liou- 
ville's theorem in that the divergence of the according 
vector field does not vanish. As we have seen, the set of 
Eqs. (J3j is not volume preserving. Consequently the so 
modified equations of motion (j3]) are no longer Hamilto- 
nian even though they produce the same identical trajec- 
tories as the Hamiltonian equations of motion (|6|) if one 
only starts with initial conditions with the same energy 
value as the one which replaces the Hamiltonian in the 
modified equations of motion. For any initial condition 



2 



with a different energy, however, the Hamiltonian equa- 
tions (J6j> and the modified equations Q will typically 
lead to different trajectories. The formal reason why the 
modified equations of motion do not satisfy Liouville's 
theorem lies in the fact that the divergence of a vector 
field at any point of phase space "sees" the full neighbor- 
hood of this phase space point and also takes into account 
directions transversal to the energy hyper-surface. 

Hoover [l[ also considers another Hamiltonian (setting 
now r = 1 and T = 1): 

# Another = («/2) [ (p / ' sf + q 2 + p* } + S ln(s) - S (7) 

and shows that for a specific initial condition yielding 
the zeroing of -ff Another, h gives the same trajectory as 
the Nose-Hoover equations. Exactly as in the case dis- 
cussed above, this is not sufficient to guarantee that the 
Nose-Hoover equations are Hamiltonian, because just like 
-HDettmann, -^Another does produce the same trajectories 
as the Nose-Hoover equations only on a specific energy 

shell ^Another = 0. 

In the last two paragraphs of his comment [l[ , Hoover 
emphasizes that the Nose-Hoover equations, when em- 
ployed to simulate heat transport, give multifractal 
strange attractors in the full system+thermostat phase 
space. Since Hamiltonian flows cannot possess attractors, 
this fact provides further evidence that the Nose-Hoover 
equations are not Hamiltonian. 

We agree with Hoover's statement, in the last para- 
graph of his comment [l[ , that if our method "is more eas- 
ily matched in laboratory experiments then it is indeed a 
step forward" . This is indeed the case. The Hamiltonian 
that we propose reads: 

ff = E^: + ^) + S + Tln X + ^ x ) (8) 

i 

where /i(q, X) is a weak interaction term between the 
system (q, p) and a log-oscillator (X, P), and b > sets 
the length scale of the log-oscillator. Its major advantage 
is that it can be matched in laboratory experiments, be- 
cause it presents the kinetic term in the standard form 
Y.iPil 2m i + P 2 /2M. In contrast, both Nose's Hamil- 
tonian ([T]), Dettman's Hamiltonian ([5]), as well as the 
Hamiltonian in (JJJ involve non-standard kinetic terms 
containing the log-oscillator position s, which can not 
readily (if at all) be matched in laboratory experiments. 
Specifically, Nose's Hamiltonian ((TJ) contains the non- 
standard kinetic term p 2 /s 2 , whereas Dettman's Hamil- 
tonian ([5|) and the Hamiltonian in ([7]) contain the non- 
standard kinetic terms p 2 / s, sv 2 /2t 2 . As we have elu- 
cidated with our manuscript |2j, one way to implement 



the Hamiltonian <j8j) is to put N neutral atoms and one 
charged ion in a box in presence of the 2D logarithmic 
Coulomb potential (this can be generated by a thin long 
charged wire). Through short ranged collisions, the ion, 
acting as thermostat, will thermalize the gas of neutral 
atoms to a Gibbs distribution of temperature T given by 
the strength of the logarithmic potential (which is pro- 
portional to the charge on the ion and to the linear charge 
density on the wire). Moreover, logarithmic potentials 
can also be created artificially by means of properly en- 
gineered laser fields for cold atoms 0] and colloidal par- 
ticles Q, and occur naturally in various situations: For 
example it is known that the motion of stars in ellipti- 
cal galaxies is governed by logarithmic potentials @ , and 
also that a logarithmic potential determines the interac- 
tion of vortices in unsteady flows ■ 

On a mathematical level, the most distinctive feature 
of our method is that, at variance with the approach de- 
scribed by Hoover in his comment [l| in our approach 
one need not choose specific values of the energy in order 
to obtain thermostated dynamics. Our Hamiltonian ([8]) 
produces thermostated dynamics globally, i.e., over the 
whole phase space: for this reason we find more appro- 
priate to qualify our thermostat as a "truly Hamiltonian 
thermostat" rather than just "Another Hamiltonian ther- 
mostat" as in Hoover's comment [l[. 

We conclude with the remark that the purpose of our 
work Q is not to criticize or to improve over the Nose- 
Hoover thermostat, whose importance and usefulness is 
out question. 
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